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Abstract 

In the numerical simulation of certain field theoretical models, the 
complex Langevin simulation has been believed to fail due to the vio- 
lation of ergodicity. We give a detailed analysis of this problem based 
on a toy model with one degree of freedom (5 = —f3 cos 9). We find 
that the failure is not due to the defect of complex Langevin simulation 
itself, but rather to the way how one treats the singularity appearing 
in the drift force. An effective algorithm is proposed by which one can 
simulate the 1/(3 behaviour of the expectation value < cos 9 > in the 
small (3 limit. 



^ Work supported in part by the Deutsche Forschungsgemeinschaft, Az. Schu95/7-2 



1 Introduction 



The complex Langevin simulation is a challenging idea for quantum systems 
with complex actions [l], @], where the conventional Monte Carlo methods are 
not efficient. It has been applied to a variety of physical systems such as QCD 
at finite temperature and chemical potential 0, ^, H, ||, gauge theories with 
external charges or fermions ||, ^, |10[. Big successes could be obtained in 
case the complex Langevin simulation works. In some cases, however, it fails 
to give correct results. 

One typical reason for the failure of the complex Langevin simulation are 
instabilities, where blow-up solutions are observed due to the non-positive 



definite Fokker-Planck Hamiltonian [11,0. To tackle this problem, the role 



of a kernel is investigated and some progress has been made |T^, |T^, [l^] . 
In principle, if a suitable kernel is introduced, blow-up solutions could be 
avoided and correct results may be obtained. Recently some more discussions 
about the convergence of the complex Langevin equation [^, |I6| and its 
spurious solutions [|l7l have been presented. 

In this paper another type of failure is concerned. In the pratical simu- 
lation of certain models, no blow-up solutions are observed but the results 
are not correct. One important example is the non-Abelian gauge theory 



coupled to external charges |18| . The violation of ergodicity due to the so- 
called segregation theorem has been blamed for this failure . This defect 
may not be remedied by the introduction of a kernel (at least up to now no 
suitable one has been found) Jl^l- It has been pointed out that the singular 
points corresponding to the zeroes in the probability distribution may play 
an important role in this problem [|19| . 

The purpose of this paper is to present a detailed analysis of the latter 
problem above. In the next section a brief review in the context of a simplified 
model with one-degree of freedom is given. Our new observation to the failure 
and some theoretical supports are discussed in section 3 and 4. In sections 5 
and 6 the possible algorithms for simulating the 5-type singular drift forces 
are presented. Finally the conclusions follow. 



2 Brief review of the problem 
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2.1 The idea and its failure 



Let us consider the following simple integral, 

^ g ^ _ /o " d9 cos^ 9 exp{/3 cos 9} _ W) 1 

/o^" d9 cos ^ exp{/5 cos 9} h{(3) (3 ' ^ ' 

The analytical result shows that the quantity is singular (~ 1//?) in the 
limit /3 — ^ 0. In this limit the denominator becomes very small, therefore, to 
simulate the quantity to high precision by calculating numerator and denom- 
inator separately becomes very difficult. In order to overcome this difficulty, 
an idea is to apply the Langevin method. For this purpose one first rewrites 
the integral as 

/■^" .r. /n^^rf^cos^expl-^effl 
< cos^^ >= / d9cos9P,ff{9) = ' 

where 

Peff oc e-^-^ff , Seff = -(3 COS 9 - In (cos 0). (3) 

This action is complex for cos 6' < 0. As the distribution P^ff oc e'^^^f is 
not positive definite, usual Monte Carlo method can not be applied. On the 
other hand the Langevin equation 

^ = -/3 sin 6^ - tan 6* + r/ (4) 

can in principle be solved and may be used for the simulation of (0). Unfor- 
tunately, this Langevin equation turns out to fail completely to generate the 
desired configurations for /5<2.5, and therefore cannot give the correct 
behaviour. This can be seen in Fig. |^. 

2.2 Conventional interpretation of the failure 
and the segregation phenomenon 

In order to give the 1/(3 singularity, the solution of the Langevin equa- 
tion 9 = 9r + i9i should get a big imaginary part. This is clear because 
^cos9 = cos 9r cosh 9i and cos 9r < 1. Contrary to this expectation, an 
important observation given in |l^ is that the configurations obtained 



through the updations by use of the Langevin equation (H) are almost real. 
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This phenomenon is always observed independently of which initial configu- 
ration one takes, and was called "the collapse of the complex distribution to 
the real distribution". 

This is surely the direct reason why the Langevin simulation fails to give 
the l//3-singularity at /? ^ as is shown in Fig. [l|. 

In addition, in |]18[ the segregation theorem has been applied to explain 



the wrong result obtained by the Langevin simulation. Due to the collapse 
of the complex distribution to a real one we can restrict the discussion to 
the stochastic process on the real axis. Let us assume that there exists 
a domain D on the real axis with a boundary on which the probability 
distribution vanishes. In this case, the segregation theorem states that the 
probability that a real diffusion process starting inside D exits the domain or 
stays permanently near the boundary of the domain, is zero. In case of our 
model (^, the zeros of Pg// divide the total configuration space V = [0,27r] 
(note that we will adopt periodic boundary conditions) into two subdomains: 
Di = [0, |)©(^, 27r] and 1^2 = (f , If the segregation theorem is applied, 
the Langevin simulation should give one of the following results, 

fn COS 6e~^^ff , , 

< cos ^ >j= 5 , i = \0T2. (5) 

depending on the starting point of the simulation. 

Note that the the above integrals have the same form as (^, while the 
integration regions are restricted to the subdomain Dx or D^- These results 
have also been plotted in Fig. |l| by dotted and dot-dashed lines. In this 
figure one may see the clear discrepancy between the result of the Langevin 
simulation and the above prediction from the segregation phenomenon. The 
situation was quite different in [0. In that paper the simulated result and 
the prediction given by the segregation assumption (especially < cos 6' >i) 
agree quite well, by which the authors concluded that the failure of their 
simulation is connected to the segregation phenomenon or, in other words, 
the violation of the ergodicity. 



3 A new observation 
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3.1 It is not the segregation but! 

Before discussing the difference between tlie result of the Langevin simulation 



obtained here and that of [|18[, let us first describe the algorithm used in more 
detail. 

In numerically solving the Langevin equation, it sometimes happens that 
the configuration comes near to the singular points where the drift force 
becomes very big. To get an accurate numerical solution in such a case, one 
should adjust the fictitious time step. Here we have followed the way used 
in ||18|. Let us write the discretized Langevin equation as A6' = D{9)At + 
v2At^ (< >= 1). The criterion to get an accurate numerical solution is 
represented by \D{e)At\ <^ |V2At^|, which gives At ^ u^e) ( have used 
C,"^ ~ 1). In discretizing the Langevin equation, we normally use a certain 
prefixed fictitious time step At = Atp. In case that the configuration comes 
near to the singular points and therefore the relation Atp < -pi^ is not 
satisfied we replace Atp by At = jjT(g^- Repeating the updations by the use 
of this adjusted time step until J2^t = '^tp, we include the last configuration 
in the set of ensemble elements. 

Let us come back to the discussion about the difference of our Langevin 
result in Fig. |I| and that in |18|. As was mentioned above the algorithm used 



here to update the configuration is almost the same as that used in |T^. The 
big difference in the result, however, arises from the number of iterations 
taken in solving the Langevin equation. In order to show this we plot in 
Fig. ^ the result of < cos 6 > versus the number of iterations. 

When one calculates the average from less than 500 iterations, the result 
coincides with < cos 6* >i, (i = 1 or i = 2) in (||), given by the segregation 
assumption. Which value (< cos^^ >i or < cos^^ >2) will result from the sim- 
ulation depends in which region, Di or D2, the starting point of the iteration 
was chosen, calculates the same quantity from a much bigger ensemble, the 
result converges to a different value that does not depend on the choice of 
the initial configuration. See Fig. |l|, which is clearly showing that the segre- 
gation is observed only within those configurations in which the number of 
updations is too small. 

This can much more directly be seen in Fig. |^, where the probability 
distribution within an ensemble of, respectively, 300, 1000 and 30000 config- 
urations are shown. In Fig. the configurations are distributed in only one 
of the regions Di and D2. On the other hand, in Fig. Ob and in Fig. lie the 
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Figure 1: 



< cos^ > simulated by the complex Langevin equation com- 
pared with the theoretical curve (solid line) given by (|I|). The 
dashed line and the dot-dashed line corresponds, respectively, to 
the theoretical prediction < cos6' >i and < cos^ >2 explained 
in subsection 2.2. The dotted line corresponds to < cos6' >abs 
explained in section 3. 



Figure 2: 

< cos 6 > simulated by the complex Langevin equation versus 
the number of iterations. Dotted line, dashed line and dot-dashed 
line corresponds respectively to < cos6' >abs, < cos 9 >i and < 
cos6' >2. 
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distribution of the configuration covers the whole region of V. 

Within the above two results, another important point needs to be stressed. 
It is the fact that the simulated results in Fig. |^ coincide exactly with|T^ 



/3 ^ f \ia QTD (Q\ /o rf^cos^exp{-5afes} 

<COSe>abs= do COS Pabs[0) = -i^— — ^ , (1) 

^0 /o dOexY>{-Sabs\ 

where 

Pabs OC 6-^""= , Sabs = -p COS ^ - lu | COS ^| . (2) 

Also note that in Fig. |c the probability distribution coincides exactly with 
Pabs, given by the line represented by circles. 

The reason why this happens will be discussed in section 4. 

3.2 What is missing then? — the (5-function 

The difference between what we should get, see Eqs.(^H), and what we have 
got, see Eqs.(|I|,^), may be found in the following simple formula 

Seff = Sabs ± i7r0(-COs6'). (3) 

where 9 is the real step-function 

From this action one gets the drift force 

=F ^7rsin6'5(cos6'), (5) 



5 Seff _ SSabs 



50 50 
and the corresponding Langevin equation 

= -psinO- P(tan 0) t sin 05{cos 0) + r^, (6) 

which is different from the normal one in Eq.(^). 

As was stated in the first paragraph of section 2.2, one will never get 
the l//3-behaviour in case the complex distribution collapses to the real one, 
independently of whether the segregation phenomenon is observed or not. 
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Figure 3: 



Probability distribution within an ensemble of configurations. 
The figure (a),(b) and (c), respectively, is that of ensembles ob- 
tained by 300, 1000 and 30000 iterations using the Langevin equa- 
tion (^). A line plotted by circles represents the distribution 
Pabs 6"'^'^''% while the solid line shows the theoretical predic- 
tion from the segregation assumption in Di. The double-circle is 
the starting point of the iterative solution. 
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An important term missing in is an imaginary part in the drift to allow 
for a complex distribution. 

In this context, the additional drift force i7i6{— cos 9) in is expected to 
work well because it kicks the configuration into the deep imaginary region. 
This pushes the distribution collapsed on the real axis to that in the whole 
complex plane and may help to give the 1/(3 singularity. 



4 Theoretical support of the idea 

In the previous section we have discussed how one may simulate correctly a 
system whose drift force is singular. Before starting to solve the new Langevin 
equation obtained there, we give some theoretical support for the idea of a 
5-function-like drift. Furthermore, an important question to be answered is 
"whether our observation contradicts to the segregation theorem or not?" . 
To make the discussion clearer we consider an alternative model which is 
simpler but has a similar structure as the cos6'-model in and d). 



4.1 A simpler model (modified Rayleigh-model) 

Let us consider the following integral, 

/+-da:x^exp(-f(x-OT 1 
^ /+~rfxa;exp(-f(x-/3)2) ^ ^ ap' 

One can rewrite it as 

J-oo jZ^dxexp{-Seff} 

where 

S.ff = ^{x-(3r-\nx. (3) 

As is seen in (|1]), the quantity < x > is also singular in the limit of /? — 0. 
This singularity arises from the fact that the denominator in (|ID vanishes in 
this limit although the numerator remains finite. The situation is very similar 
to that of the cos6'-model. Therefore, this model will be a good candidate 
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to check all of the points discussed in the previous sections. Exactly in the 
same way as in case of the cos 6'-model, the naive Langevin equation^ 

X = -a{x - P) + - + r]{t). (4) 

X 

fails to reproduce the l/j3 singularity given in (|I|). Based on this simple 
model, therefore, we will first study the segregation theorem in more detail 
in the next subsection. 



4.2 The segregation phenomenon 

Let us start from the Fokker-Plank equation equivalent to the above Langevin 
equation, 

Pix;t) = ^^[B{x,t)P{x,t)]- ^[Aix,t)Pix,t)], 

B{x,t) = 2, A{x,t) = A{x) = -a{x - p) + (5) 

This is a second-order parabolic differential equation. In order to specify its 
solution one needs an initial condition, e.g., P{x; to) = S{x — Xo), and bound- 
ary conditions. In case the drift force A{x,t) is singular (or the diffusion 
coefficient B[x,t) vanishes) at some point x = Xb, a. certain boundary con- 
dition should be imposed at this point which is called a prescribed boundary 
condition. What kind of boundary condition one should take at x = Xb is 
determined by the form of the drift and the diffusion coefficient. An exten- 



sive analysis of this problem has been given by Feller |21[]. See also ||2^ for a 
review. 

For our Langevin equation (Q) the prescribed boundary condition will 
be given at the singular point of the drift, x = Xb = 0. [xb,Xo]. Following 



the classification of Feller the singular point x = Xb = should be 



treated as an entrance boundary. In this case the solution of the Langevin 
equation (^) should be such that if the configuration starting from, e.g., 
X = xo > Xb comes near to that boundary it should be pushed back to the 



Note that this Langevin equation is closely related to the well known Rayleigh process, 
-ax+ +r){t). 
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original direction x > Xb- Therefore, the configuration starting from x > Xb 
can never leave that region^. This is the so called segregation phenomenon. 



4.3 The reason for the failure of the naive Langevin 
equationand the role of the ^-function-like drift 

Similar to what we saw in Fig. 2 and Fig. 3, however, this is not the case when 
one naively does a numerical simulation for solving the Langevin equation (^) 
by discretizing the time t. The configuration turned out to pass through the 
barrier at x = Xb easily and the segregation theorem looks to be violated. The 
main reason why this happens may come from the fact that by discretizing 
the time we can not treat the singularity in the drift rigorously. The adjusted 
time step has been introduced in order to take into account the effect of the 
singular region as rigorous as possible. But this seems still not enough. 
The system which are simulated by the discretized version of the Langevin 
equation (|) can not exactly be the process given by (^) but that with some 
regularized drift A{x)r. In Figure 4, A{x)r is taken to be constant ±l/ro for 
|x| < tq. The form of A{x)r in \x\ < tq is, by no means, unique. However, 
one can easily check that the following discussion does not depend on the 
way how one regularizes A{x) in this region. 

For the stochastic process with the above regularized drift A{x)fi, no 
prescribed boundary condition exists by definition. The solution, therefore, 
can pass through the point x = without any problem. In addition to 
this, an interesting point is that for this regularized process the stationary 
solution for P{x) (at |a;| > tq) is given by Pabs oc |x|exp(— |(a; — and 
not by xexp(— |(a; — To show this, let us try to solve the stationary 

solution P{x) in the case of the regularized drift. 



The solution can be represented as (for simplicity we take /3 = in the 
following) 



^ When the configuration starts exactly from that point x = Xb it should enter either 
the region x > Xb ot x > Xb, which is the origin of the name entrance boundary. 





To < X 



P{x) = Cie 
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< X < ro 
-ro < X < 

X < —Tq 



, 



P{x) = C3e-^^"-^^ (7) 

with ci,C2,C3 and C4 being constants. The continuity conditions of P{x) at 
X = ±ro and x = gives an interesting relation between ci and C4, namely 

C4 = cie^°('^) = -ci. (8) 

This proves that, outside of the singular region, the Fokker-Planck distribu- 
tion has the form of P{x) oc \x\ exp{(— a/2)x^}. From those discussions, it 
may be clear why our numerical simulation gave the results corresponding 
to Pahs{0) given in (|l|) and in (|) in the previous section. 

Our proposal for the way out of this problem mentioned in the last sub- 
section is to include an extra 5-function type drift in the Langevin equa- 
tion due to the relation Inx = ln|x| ± Z7r6(— x). The drift coming from 
ln|x|, i.e., — (5(ln |x|)/5x = P(l/x), is corresponding in some sense to the 
regularized drift A[x)r discussed above. Then the total drift should beQ. 
A{x) = A{x)ji + i7r5(x) It is interesting to repeat the same discussion as 
above using the new drift A{x). For this drift A{x), the solution of J = 
{dx + A{x))P{x) = can not, in general, be continuous at x = 0. This is 
seen by integrating both sides of this equation in a small region [— e, +e], and 
taking the limit e — > 0. On the other hand applying the same discussion 
to the stationary F-P equation, P = dxJ = dx{dx + A{x))P = 0, one finds 
that dxP{x) can still be continuous at x = 0. From this continuity condi- 
tion of dxP{x) one gets C4 = ci rather than that in (H). This means that 
P(x) oc xexp(— ax^/2) rather than P(x) = Pabs{,x) oc |x| exp— ax^/2. 

This shows that if we succeed to estimate the effect of the 5-function type 
drift in the Langevin equation, we will get the correct result. 

5 How can one do the simulation for the 5- 
function type drift? 



^ Note that an addition of this J-function type drift will never spoil the regularity of 
the stochastic process at a; 0. This can easily be seen by calculating the functions 
f{x), g{x), hi{x) and h2{x) defined in ||2^ by using this drift. All of them turn out to be 
integrable in the region of [0, xo] for any finite xq. 
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5.1 A model with a real 5 function 



In order to study how to solve the Langevin equations with the 5-function 
type drift numerically, we first consider an even simpler model given by the 
real action 

S2 = ^{x-f3r--e{x), (1) 

where c > is a real constant, ue to the existence of the —ce{x)/2 in the 
action, the distribution P2 oc e""^^ has a discontinuity at x = 0, see Fig. 5. 
The probability distribution is shifted to the direction of positive x, which 
increases the value of, e.g., < x >. 

In the Langevin simulation, this effect will be given by the drift —6/6x{—c/2e 
c6{x) in the corresponding Langevin equation 

X = —a{x — P) + c6{x) + r] (2) 

When the configuration crosses x = while updating the Langevin equation, 
it gets a kick towards positive x, which works to increase the expectation 
value < a; >. 

5.2 A success of the numerical simulation 

For the numerical solution of the Langevin equation we first have to fix how 
to treat the drift c6{x). In this subsection we like to discuss two ways, 

— the use of a smeared 5-function, and 

— an integration of the (5-function around x = 0. 

The use of the smeared (5-function: One can replace the 5-function 

by the smeared one, 

S{x) — , S, = - ^ , (3) 

TC x^ + 

where e is a small number. In Fig.6, the results for < x > and < x^ > 
simulated by the Langevin equation with 6{x) being replaced by the 
smeared one S^{x) are given. One can see that, within this real model, the 
effect of the discontinuity in the probability distribution is nicely simulated 
by the effect of the 5-function type drift in the Langevin equation. 



12 



Figure 4: 

(a) Regularized drift force 74(x)/j at x ~ 0. (b) Solution of (H) 
under an assumption of continuous P{x). (c) Solution of (^) under 
an assumption of continuous dxP{x). 



Figure 5: 

Schematic figure of the probability distribution P2 oc e^"^^ with 
5*2 given in (0). 
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Integrating out the (5-function: It is also possible to get an effective 
algorithm by analytically integrating out the (5-function. Let us rewrite the 
Langevin equation in a discretized form, and consider the step from Xi 
to Xf (Ax = Xf — Xi), 

Ax = a + c , At, 
Ax 

+ At, (4) 

where ^ is a Gaussian noise with < >= 1. When x does not cross x = 
within this step, AQ{x) = and therefore from (H) we simply get 

Ax = a (5) 

When X crosses x = 0, A0(x) = ±1 ( depending on whether Xj > or Xj < 
) and the Langevin equation (^) becomes a second order algebraic equation 
for Ax. Solving this equation we get 

Ax = (xi<0,a>0). 



a — \/ — 4cAt , „ . , 

Ax = (xi>0,a<0). (6) 

The simulated results for < x > by the use of this algorithm is shown in 
Fig. ^ One can see that the Langevin simulation with the 5-function type 
drift is being able to recover the exact one very nicely. 



6 Complex (5-function 

In this section we come back to our original problem, the cos^-model. As 
discussed in section 3.2 the Langevin equation for this system is given by 
(^. This Langevin equation has an additional drift i7r(5(cos^) that kicks the 
configuration Q into the deep imaginary region. This will surely help to give 
a big < cos Q > because of the reasons explained in section 2.2. The problem, 
however, is how to continue this (5-function into the whole complex ^-plane. 

In order to explain the situation more clearly, let us once more use the 
modified Rayleigh model given in (|I|). The Langevin equation for this system 



14 



Figure 6: 

The result of < a; > and < > (for a = 0.2, (3 = 0) simulated 
by the above Langevin algorithm versus c. For each run 50000 
iterations have been done with At = 0.04. Errors are estimated 
from 32 runs of this kind. The sohd line represents the exact result. 



Figure 7: 

The result of < a; > and < > (for a = 0.2, /3 = 0) simulated 
by the above Langevin algorithm versus c. For each run 50000 
iterations have been done with At = 0.01. Errors are estimated 
from 32 runs of this kind. The sohd line represents the exact result. 
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is given by with the additional drift i7TS{x). In this case, even when one 
starts the iteration from some point on the real axis, it becomes complex 
after getting a kick by the (5-function. Then the question arises how to treat 
the inSlx) in the complex plane z = x + iy. 

Smeared 5-function: One possibility may be to use the smeared 

5-function continued into the complex plane, i.e., S^{z) = e/7i{z^ + e^). But 
this does not work. The reason may be that the above smeared 6 function 
which has been introduced in order to take into account the effect of the 
singularity in the drift 1/x has again a new singularity at z = ±ie in the 
complex plane. Actually we have tried to use it for simulation, it could not 
reproduce a correct 1/(3 behaviour for < cos^ >. 

Integrating out the 5-function: We also have a problem in ap- 

plying this idea to the complex case, because we do not have the formula 
corresponding to (|^) in the complex case. But neglecting this fundamental 
question, we may use the solution of this equation, i.e., (|^) and (|^), with 
c = in. The result is not good enough to produce the singularity al- 
though it refines the results obtained without the 5-function-type drift. 

In spite of these failure, it is still very encouraging that with the idea 
of the (5-function type kick we can find an effective algorithm by which the 
above 1/P behaviour in the cos^ model is successfully simulated. Let us 
approximate the drift in sm6S{cos6) in the Langevin equation (|^) in the 



following way |19]: Consider the region T) in the complex plane of ^ = 9r+i9i 



around the singular points 6 = it/ 2 or 37r/2, 

6 eV ^ \^cose\ < S, (1) 
6 being a small constant, and take 

In Fig. 1^ we show the result of < cos 9 > for the cos 9 model. The data drawn 
by black squares are those simulated by the use of the Langevin equation (|^), 
where the drift proportional to (5(cos 9) has been replaced by the kick given 
by equation (0) above. One can see that the simulated result coincides with 
the theoretical prediction quite nicely in the whole region of jS. 

In spite of this wonderful success of the numerical simulation, we have not 
yet found any smart theoretical justification of the above algorithm. More- 
over the effectiveness of this algorithm can not be universal for other models. 
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Figure 8: 



The black square and the circle represents, respectively, the re- 
sult of < cos 9 > with and without the kick (^. For the parameters 
6 = 0.07,At = 0.002 have been taken. In order to calculate the 
average 100000 iterations has been done for each run and the error 
has been estimated from 8 different kind of runs. The solid line 
represents the theoretical prediction. Dashed line and dotted line 
corresponds, respectively, to the theoretical prediction < cos 6* >i 
and < cos^ >abs explained in subsection 2.2 and section 3. 
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The modified Rayleigh model discussed in section 4 is a good example. If the 
5-function type drift in6{x) for this model is treated exactly in the same way 
as above, it will just give the imaginary kick iir to the configurations near 
to the singular point z ~ 0. The big imaginary part in z, however, does not 
directly give a big ^ < z >=< x > contrary to the case of the cos^ model. 
We need therefore deeper understanding of the algorithm found above from 
more theoretical view point like that in the case of those models with real 
5-function in the previous section. 

7 Conclusion 

In the numerical simulation of certain field theoretical models, the complex 
Langevin simulation has been believed to fail due to the violation of the er- 
godicity. In this paper we have given a detailed analysis of this problem based 
on a toy model in one degree of freedom (cos 9 model) . The corresponding 
Langevin equations involved in the above problem have a singular drifts, e.g., 
oc tan^. Our observation is that the failure is not due to the defect of the 
complex Langevin simulation itself, but rather to he way how one treats the 
above singular drift force. We have tried to justify this statement using also 
some alternative models. Under the above observation we could also give an 
effective algorithm by which we can simulate wonderfully the behaviour 
of the expectation value < cos^ > in the limit of /5 — > 0. 

Unfortunately, however, the final theoretical justification of the above 
effective algorithm is still missing. When one succeeds in getting the rigorous 
theoretical background about how to treat the singular drift in the complex 
Langevin simulation, practical gain in the numerical simulation of lattice 
field theory will be extremely large. 
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